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ABSTRACT 

We estimate minimum dust abundances required for secular gravitational instability (SGI) to op- 
erate at the midplane dust layer of protoplanetary disks. For SGI to be a viable process, the growth 
time of the instability T grow must be shorter than the radial drift time of the dust T drift . The growth 
time depends on the turbulent diffusion parameter a, because the modes with short wavelengths are 
stabilized by turbulent diffusion. Assuming that turbulence is excited via the Kelvin-Helmholtz or 
streaming instabilities in the dust layer, and that its strength is controlled by the energy supply rate 
from dust accretion, we estimate the diffusion parameter and the growth time of the instability. The 
condition T grow < Thrift requires that the dust abundance must be greater than a critical abundance 
Z m [ n , which is a function of the Toomre parameter Q g and aspect ratio h g /r of the gas disk. For a 
wide range of parameter space, the required dust abundance is less than 0.1. A slight increase in dust 
abundance opens a possible route for the dust to directly collapse to planetesimals. 
Subject headings: planets and satellites: formation — protoplanetary disks 



1. INTRODUCTION 

Formation processes of planetesimals have not been 
well understood in planet formation theory. The grav- 
itational instability of the dust layer at the midplane 
of a protoplanetary disk has been proposed as a pos- 
sible route to planetesimal formation. The classical sce- 
nario of gravitational instability requires high densities 
of dust layers to surpass the Roche limit (Goldreich & 
Ward 1973; Sekiya 1983). Dust sedimentation only in 
the vertical direction hardly achieves such a high density 
(Sekiya 1998). The dust layer becomes turbulent via 
Kelvin-Helmholtz (KH) or streaming instabilities when 
its dust density exceeds the gas density, which is much 
less than the Roche density, and further dust settling is 
significantly suppressed (Chiang & Youdin 2010 for re- 
view). Radial drift of the dust caused by gas drag may 
provide a possible route for further accumulation of the 
dust. When the dissipative effects of gas drag are in- 
cluded, the dust layer is secularly gravitationally unsta- 
ble to the modes of dust accumulation in the radial di- 
rection (Ward 2000; Youdin 2005a, 2005b). This secular 
gravitational instability (SGI) occurs in any dust layers 
even if their densities are small (i.e., there is no crite- 
rion on their Toomre's Q values for SGI; Youdin 2011; 
Shariff & Cuzzi 2011; Michikoshi et al. 2012; these pa- 
pers are referred hereafter as Yll, SC11, and MKI12, 
respectively) . 

However, the growth rate of SGI is greatly suppressed 
by gas drag, particularly for small particles, and could 
be too slow to have any effect on planetesimal formation. 
Yll, SC11, and MKI12 argue that the growth timescale 
must be shorter than the lifetime of the dust drifting 
toward the star owing to gas drag. The wavelength of 
the most unstable mode is determined by the balance 
between the radial accumulation of dust particles due 
to self-gravity and their diffusion due to gas turbulence. 
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Only perturbations with wavelengths sufficiently long are 
unstable. For stronger turbulence, wavelengths of the 
unstable modes, as well as growth timescales, are longer. 
Yll and SC11 have calculated the growth timescales for 
various values of the turbulent diffusion parameter a, 
and have derived an upper limit on a required for the 
growth timescale to be shorter than the drift timescale. 
Their approach is general and can be applied for any disk 
turbulence, provided that the value of a is known. 

In this paper, we focus on disk turbulence induced in 
the dust layer, assuming that the gas disk itself is a glob- 
ally laminar flow. In such a disk, the concentration of 
dust particles at its midplane triggers turbulence. The 
velocity difference between the dust and gas induces the 
KH or streaming instabilities and causes turbulence in 
the dust layer (Chiang & Youdin 2010). The ultimate 
energy source for turbulence is the accretion energy of 
the dust drifting toward the star, regardless of the kind 
of instability that occurs in the dust layer. Thus, a 
simple energetics can be applied to estimate the turbu- 
lence strength, provided that the dust accretion velocity 
is known. Takeuchi et al. (2012, hereafter T12) calculate 
the energy supply rate to turbulence, using the classical 
formulae on the drift velocity of the dust (Nakagawa et 
al. 1986; Weidenschilling 2003; Youdin & Chiang 2004), 
and then estimate the turbulence strength or the value 
of a. T12 has shown that, from the comparison between 
the estimated value of a and the recent results of numer- 
ical simulations of turbulence by Johansen et al. (2006) 
and Bai & Stone (2010), the turbulence strength excited 
via KH and/or streaming instabilities can be estimated 
from the dust accretion rate. 

While Yll, SC11, and MKI12 have derived the growth 
time of SGI as a function of a and other disk parameters 
(Equation (51) of Yll), T12 has obtained a values of tur- 
bulence in the dust layer (Equation ([5]) below) . Combin- 
ing these results gives an expression of the growth time 
without the unknown parameter a, as shown in Section 
12.11 Consequently, in Section [231 the criterion for SGI 
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to operate is obtained as a condition on the dust abun- 
dance Z in the disk. For SGI to operate faster than the 
dust drift timescale, the dust abundance must be greater 
than a specific critical value Z m in, which is a function 
of the Toomre Q value of the gas disk and the deviation 
fraction r\ of the gas velocity from the Keplerian velocity. 
The Z m i n value derived in this paper gives the minimum 
dust abundance required for SGI, because we consider 
only turbulence induced in the dust layer. If turbulence 
were caused by another mechanism such as magneto- 
rotational instability (MRI) of the gas disk, a would be 
larger than the value we adopted, and consequently the 
required value of Z m - m would increase. Hence, Z m - m de- 
rived in this paper is considered as the minimum value 
required for SGI in realistic disks. 

2. CONDITION FOR SECULAR GRAVITATIONAL 
INSTABILITY 

We consider a protoplanetary disk initially in a laminar 
flow state. Sedimentation of dust particles to the mid- 
plane of the disk induces turbulence via KH or streaming 
instabilities caused by velocity differences between the 
dust and gas. A (quasi-)steady dust layer forms when 
turbulent diffusion matches dust settling, in which steady 
state turbulence is maintained. We neglect intermittency 
of turbulence. If an extra source of turbulence is present 
such as global turbulence of the gas disk via MRI, the 
turbulence would be stronger than that excited only by 
the dust-gas velocity difference. Thus, in this paper, we 
consider the minimum strength of turbulence. 

The dust particles are characterized by their stopping 
time t s , which is the timescale of damping the velocity 
difference from the gas. We define the non-dimensional 
stopping time, T s = t s flK, normalized by the Keplerian 
frequency £!k- In this study, we focus on the dynamics of 
relatively small particles such that T s < 1. Particles with 
T s ~ 1 experience high speed collisions and the fastest 
radial drift, which likely hinder particle growth. SGI of 
the dust layer is a possible route for directly forming large 
bodies (T s 1) from small particles (T s 

In the following subsections, we compare the growth 
time of SGI with the orbital drift time of the dust parti- 
cles due to gas drag in order to determine the condition 
for SGI to be a relevant process for planetesimal forma- 
tion. 

2.1. Timescale of Secular Gravitational Instability 

In his Equation (51), Yll shows that the growth time 
of SGI for small particles with T s < 1, normalized by the 
Keplerian time ^l^ 1 , is 
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Here a is the turbulent diffusion parameter (see below), 
Q g is the Toomre stability parameter of the gas disk, 
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growth time depends on the turbulent diffusion parame- 
ter a through its stabilizing effect for perturbations with 
short wavelengths. 

Turbulence in the dust layer is induced by velocity dif- 
ferences between the dust and gas. If the drag force be- 
tween the gas and dust were not effective, the dust would 
orbit with the Keplerian velocity «k, while the gas would 
orbit with a sub-Keplerian velocity 

v 9 = C 1 - v)vk , (3) 

where 77 = -(2p g rQ^)- 1 dP '/ 'dr ~ (c g /v K ) 2 ~ 1CT 3 - 
10 ~ 2 , p g is the gas density, and P is the gas pressure. 
In the calculation of a, a slightly different definition of 
77 is useful (see T12 or Appendix |A1 in this paper). We 
introduce 77 defined by 
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V 2 = CvV , 
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where C v is a factor of the order of unity, which depends 
on the radial profile of the gas disk. In this paper, we 
adopt C n = 1 (i.e., fj = rj) for simplicity, except in Section 
[3] where specific disk models are considered. 

T12 shows that the dust particles accrete toward the 
star due to gas drag and supply energy to turbulence 
via KH and/or streaming instabilities. Dust accretion is 
either caused by the gas drag acting on individual par- 
ticles or by turbulent drag acting on the surface of the 
dust layer. The former (individual drag) is effective if the 
dust-to-gas ratio at the disk midplane, / m i<j = Pd/Pg\ z=0 , 
is less than unity, and the latter (collective drag) is ef- 
fective if /mid ^ 1- Appendix [A] briefly summarizes dust 
accretion due to both drag types (see T12 for detailed 
discussions). The turbulence strength or the parameter 
a is determined by the energy supply rate from dust ac- 
cretion toward the star. The approximate expression of 
a for T s < 1 particles is given by 



(5) 



where C e e = 0.19 is the energy supply efficiency (see 
T12), and C\ = 1.0 and C2 = 1.6 are the numerical 
factors. This expression connects the approximate for- 
mulae of a for the two limiting regimes of / m id <C 1 and 
/mid > 1 (Equations dA"2|) and (|A"3j) ). Note that the 
condition / m id <C 1 (or / m id 3> 1) corresponds to the 
condition Z < (C cff ?7) 1/2 (or Z > (C c ff?7) 1/2 ) (see Equa- 
tion (|A5p ). The numerical factors C\ and C2 are ad- 
justed to make an appropriate fit to the overall behavior 
of the numerical result. A comparison of this approxi- 
mate expression to the numerically calculated a is shown 
in Appendix [A] Substituting Equation ([5]) in Equation 
(HI) gives 



Q 
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Z 2 T S 



(dC^fjZ) S + (CaCeff^Z- 1 ) 



where c g is the sound speed of the gas, G is the gravita- 
tional constant, and S g is the surface density of the gas 
disk. The disk "metallicity" Z is the ratio of the surface 
densities between the dust and gas, Z = S^/Sg. The radial drift of dust particles due to gas drag is estimated 



( 6 ) 

showing that T grow is a function of the disk parameters 
(fj, Z, Q g ) and is inversely proportional to T s . 

2.2. Timescale of Radial Drift 
According to T12, the non-dimensional timescale of the 
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(a) log(Z m/n ) for 

Tgrow < 7 drift 




Fig. 1. — Ratio of the growth time to the radial drift time, 
Tgrow/Tdrif t , against dust abundance Z. (a) For models with a 
fixed Toomre parameter of the gas disk, Q g = 30. (ft) For models 
with a fixed deviation fraction of the gas velocity from the Keple- 
rian velocity, r\ = 3 X 10 — 3 . The horizontal dotted line indicates 

Tgrow — Thrift- 
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where the midplane dust-to-gas ratio / m id is approxi- 
mately given by 



/mid 



CiC ctf rj 



z 2 



C 2 C c fffi 
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In the above estimate, both individual and collective 
drags are considered. Derivation of the above expres- 
sions are described in Appendix [A"l The timescale of the 
radial drift is a function of (77, Z) and is inversely pro- 
portional to T s . 

2.3. Comparison of Timescales 
2.3.1. Minimum Dust Abundances 

Yll, SC11, and MKI12 show that dust layers are al- 
ways unstable and subject to SGI, but the growth rate for 
particles tightly coupled to the gas is strongly suppressed 
by gas drag. The growth timescale can be larger than the 
other timescales, in which the dust layer evolves signif- 
icantly by other processes such as the dispersal of the 
gas disk and the radial drift of the dust to the star. The 
condition for SGI to be a relevant process for planetesi- 
mal formation dictates that its growth timescale must be 
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Fig. 2. — Contour map of the minimum dust abundance Z m i n on 
the fj-Qg plane, (a) Z m ; n required for the condition T grow < 7d r ift. 
The contours are labeled by log(Z m ; n ). The dashed line shows the 
locus of Z min = (Ccff^) 1 / 2 (i.e., / mid fa 1 when Z = Z min ). (6) 
Z m i n required for the condition 10T g row < T^ft. 

shorter than the other evolution timescales. Yll shows 
that the radial drift imposes the most stringent condition 
for a wide range of disk parameters. We discuss the con- 
dition for T grow < Tdrift- Because both T grow and Tdrift 
are inversely proportional to T s , the ratio T gr0 w/Tdrift is 
independent of T s ; that is, the condition Tg r0 w < Tdrift is 
not affected by particle size. We consider dependence of 
Tgraw Adrift on ^ e parameters (77, Z, Q g ). 

First, we examine the limiting cases of Z <C (C ff??) 1//2 
and Z > (C e ff?7) 1/2 (/mid < 1 and / mid > 1). In both 
limits, 



T 

J- cri 



Tdrift 



(dCoff ) 2/3 rf/^z-Va for z <K (Ceefj) 1 / 2 
2C- 1 (C 2 C cS ff, 4 Q 2 g Z- 6 for Z » (C^r?) 1 / 2 



(9) 

For Z < (Ccffrj) 1 / 2 , Tgrow/Tirift is proportional to 
Z~ 4 /' 3 , while for Z » (Cegfj) 1 ^ 2 it rapidly decreases as 
Tgrow /Tirift cx Z~ 6 . In Figure □] T g row/T dr ift is plot- 
ted against Z for various values of 77 and Q g . For 
Z > (CoffT)) 1 / 2 ~ lO- 2 , it is evident that T &ovl /T dm 
rapidly decreases. For sufficiently large values of Z, 
Tgrow/Tir ift is less than unity, indicating that SGI oper- 
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Fig. 3. — Contour map of the growth time multiplied by the stop- 
ping time, T grow T s , on the fj-Q g plane. The contours are labeled 
by log(T grow T s ). In calculating T grow (»?, Q g ), the minimum dust 
abundances Z m i n (fj, Q g ) are used. The dashed line shows the locus 

of Z min = {C eS fi) 1/2 - 

ates faster than the radial drift of the dust. Prom Figure 
[TJ the minimum dust abundance required for SGI, Z mYn , 
is determined as Z satisfying Tgrow/Tdrift = 1. The min- 
imum dust abundance Z m in ranges between 10~ 2 and 
10 _1 for parameters (fj,Q g ) adopted in Figure [TJ 

In Figure [Jf a), the minimum dust abundance Z m i n re- 
quired for Tg r ow < Tdrift is plotted as contours on the fj- 
Qg plane. The locus of Z min = (Cegfj) 1 ^ 2 is represented 
as a dashed line on the fj-Q g plane. On this line, the 
midplane dust-to-gas ratio / m id becomes close to unity 
when Z = Z m in- Substitution of Z = (Cesfj) 1 / 2 into 
Equation ([9]) shows that the condition Tgrow/Tdrift = 1 



(i.e., Z = Z m in) becomes Q g 



V 



- 1 / 2 . Thus, the dashed 



line in Figure [J^a) is linear with a slope of —1/2. Above 
this line, Z min > (C cS fj) 1/2 (and / mid > 1 for Z = Z min ); 
below it, Z m in < (Ccff^) 1 / 2 (and / m id < 1). The depen- 
dence of Z m - m on (fj, Q g ) is expressed as 



(SC^ClC^fQpV* for 
{IC-^ClC^Ql,)^ for 



h « r 1/2 
} a » r 1/a 



. ( 10 ) 

Larger dust abundances are required for larger fj (i.e., 
hotter gas disks) and larger Q g (i.e., less massive gas 
disks). In standard models for protoplanetary disks, fj = 
10~ 3 - 10^2, and Q g = 10- 1000 (see Figure © below). 
Thus, if the dust abundance Z is larger than 0.1, such 
dusty disks operate SGI. If the gas disk is so cold that 




Fig. 4. — Contour map of the diffusion parameter divided by the 
stopping time, a/T s , on the fj-Qg plane for the dust layer with 
Zmin- The contours are labeled by log(a/T s ). The dashed line 
shows the locus of Z mm = (Ccfffj) 1 / 2 . 

disks (Qg ^> (i???) -1 / 2 ; above the dashed line in Fig- 
ure [2]), Zmin depends only weakly on R. For massive 
gas disks (Q g -C (Tfty) -1 / 2 ), required Z mm is higher for 
larger R (Z mm oc R 3 ^ 4 ), but would not be too high (i.e., 
-^min ;$ 0T even for R = 10; see the region below the 
dashed line in Figure HJb)). In the following discussions, 
we use the condition T grow < Thrift- 

In this paper, it is assumed that about 20% of the 
dust accretion energy is used for turbulence excitation 
(Ccff = 0.19). This value is determined by compari- 
son with the numerical simulations of turbulence excited 
via KH instability by Johansen et al. (2006). Because 
their simulation was two-dimensional, the realistic value 
of C e ft could be different. For example, three-dimensional 
simulations by Lee et al. (2010) show that the critical 
Richardson number for KH instability increases with Z , 
implying that C c g also increases with Z (see discussion 
in Section 5.2.1 of T12). As shown in Equation fllT " 



1/2 

C^r. For the maximum efficiency (C, 



Zmir 

Z m in would be about 2 times larger than our estimate. 



off 



1), 



2.4. Conditions on T gr 



grow ; 



and A 



In the previous subsection, the condition for Tg row < 
Tdrift was discussed. Next, we discuss other constraints 
required for SGI operation. First, we consider the con- 
dition for the growth time to be less than the disk life 
time. The growth time multiplied by the stopping time, 
Tg m wTs, is plotted on the fj-Q g plane in Figure [3] The 



fj is as small as 10 and if it is so massive that Q g is values of T gIOW (fj, Q g ) are calculated using the minimum 
as small as 10, then even a standard value Z = 0.01 is 
sufficiently large to operate SGI. 



2.3.2. Dependence of Zmin on Some Model Parameters 

The growth time T grow is just a time for density per- 
turbations to increase e-fold. Because many e-folds are 
needed for a significant density increase, actual condi- 



tion for SGI would be RT S , 



< Tdrift, where R > 



1 is the required number of e-folding. From Equa- 
tion (j9]) the minimum dust abundance Z m in scales as 
Z min oc i? 3 / 4 for Q g < (i?7 ? )" 1 /2 and Z mln oc R 1 ^ 
for Q g » (Rrj)' 1 / 2 . Figure [2fb) shows Z min derived 
for the condition 10T grow < Tdrift. For less massive gas 



dust abundance Z — Z lrL i n (fi,Qg), the values of which 
differ according to (fj, Q g ) as shown in Figure [21 For 
larger values of Z, T grow is shorter (T gr0 w oc Z~ 4 / 3 for 
Z < (C cff fj) 1 / 2 and T grow oc Z~ 4 for Z > (C cS fj) 1/2 ). Be- 
cause T g row is inversely proportional to T s , the contours 
are labeled by log (TgrowTs). Thus, the growth timescale 
is obtained by dividing the value in Figure [3] by T s . For 
fj > 10~ 3 and Qg < 10 2 , T grow is less than lO 3 ^ 1 ^ 1 . 
Thus at 1AU, the dust layer composed of T s > 10~ 3 
particles (size a > 1mm) with the dust abundance Z m i n 
operates SGI within a disk life time of ~ lMyr. 

Figure [4] shows the diffusion parameter a due to tur- 
bulence induced in the dust layer with the abundance 
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Zmin- Because a is proportional to T s , the contours are 
labeled by log(a/T s ). Thus, the a value is obtained by 
multiplying the value in Figure |4] by T s . For ij > 10~ 3 
and Q g < 10 2 , a is larger than 10 _4 T S , except for very 
small Q g < 10 and ij < 3 x 10~ 3 . If the disk is turbu- 
lent due to other mechanisms than that originating from 
the dust layer and if turbulent diffusion is stronger than 
that shown in Figure El then the required value for Z 
would be higher. Even in the dead zone where MRI is 
not active, the gas could have turbulent motion, causing 
diffusion of the dust (Fleming & Stone 2003). Okuzumi 
& Hirosc (2011) show that lower values of the diffusion 
coefficient in the dead zone are realized for a wider dead 
zone or weaker vertical magnetic fields. For example, if 
the plasma j3 (the ratio of gas pressure to magnetic pres- 
sure) is greater than 3 x 10 6 , the diffusion coefficient due 
to MRI turbulence is less than 10~ 5 (see model Xlb of 
Okuzumi & Hirose 2011). In such a weakly magnetized 
disk and for T s > 0.1 particles, turbulence originating 
in the dust layer is stronger than MRI turbulence in the 
dead zone. 

The wavelength of the most unstable mode must be 
smaller than the disk radius. The wavelength A is given 
in Equation (56) of Yll. The ratio of A to half of the 
disk radius r is 

r C 7 jZT s 

and is plotted in Figure [5] for disks with the minimum 
dust abundances Z m i n . Note that 2A/r does not depend 
on T g , because a is proportional to T s . Figure [5] shows 
that A is smaller than r/2 if 77 is less than 3 x 10~ 2 , al- 
though it is still comparable to r/2 even for fj as small 
as 3 x 10~ 3 . In the analysis by Yll, SC11, and MKI12, 
the gas disk is assumed to behave as a stationary back- 
ground, and its velocity profiles are not affected by the 
gravity of the dust layer. However, the response of the 
gas to perturbations longer than the gas scale height is 
not clear. For perturbations with smaller wavelengths to 
be unstable, larger values of Z than those shown in Fig- 
ure [2] are required. The wavelength varies as A cx Z~ 1//3 
for Z < (Ccff?]) 1 / 2 and A cx Z~ 3 for Z > (C eS fj) l/2 ■ In 
the latter regime (Z > (Cogfj) 1 ^ 2 ), an increase by a factor 
of 2 in Z results in an order of magnitude decrease in A, 
while in the former regime, reduction in the wavelength 
requires a large increase in Z. 

3. MINIMUM DUST ABUNDANCES FOR VARIOUS DISK 
MODELS 

In the previous section, minimum dust abundances 
for SGI were obtained for given parameters (fj,Q g ). In 
this section, we consider several disk models. We adopt 
power-law profiles for density and temperature distribu- 
tions of the gas disk, such that 

S 9 = S 9,0?AU > ( 12 ) 

T = Tor-Xxj , (13) 

where tatj is the distance from the star measured in AU. 
In such a disk, Q g cx rj^ -9-3 ^ 2 and fj cx r 1 ^ . We 
consider two cases of hot and cold gas disks, To = 300K 
and 150K, while we fix the power-law index q = 1/2. 
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Fig. 5. — Contour map of the normalized wavelength of the most 
unstable mode 2A/r on the fj-Q g plane for the dust layer with Z mnl . 
The contours are labeled by log(2A/r). The dashed line shows the 
locus of Z min = (Ceffr/) 1 / 2 . 
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Fig. 6. — Models of the gas disks are plotted on the fj-Q g plane. 
Red and blue lines represent models with shallow and steep density 
profiles, E oc r _1 and S <x r -1,5 , respectively. Solid and dashed 
lines indicate models of cold and hot gas disks, To = 150K and 
To = 300K, respectively. Dots on the lines show locations on the 
disks, 0.1, 1, 10, and 100AU. The disk mass is M disk = 0.03A%. 
Contours show log(Z m j n ). 

For the density profile, p = 1.0 and 1.5 are adopted and 
the mass of the gas disk inside 100 AU is varied within 
Miisk = 10~ 2 — 1O _1 M . Figure [6] shows the variability 
in the values of fj and Q g between the models and with 
r. In this figure, models with Mdisk = 3 x 10~ 2 M Q are 
shown. For various values of Mdisk, Qg scales as Q g cx 
M disk . 

Figure [7] shows the required minimum dust abun- 
dances Z m { n against tau for various disk models. For 
the steep density profile (p = 1.5; blue dashed lines), 

Zmin cx r A {^. (In the regime of Q g <C f]^ 1 ! 2 of Equa- 
tion Q3JJ, Z min cx ?y 5 / 4 Qg /2 cx r A 7 u; in the Q g > fj- 1 / 2 

regime, Z m - m cx ff^ 3 Qg cx rly). SGI is more viable 
at the inner part of the disk. For p = 1.0 (red solid 

lines), Z min cx r A y 2 only weakly depends on r (in the 
Q g r/~ 1//2 regime). The outer part of the disk with 
Mdisk = 1O _1 M and p = -1 is in the Q g < i)- 1 ' 2 
regime, and Z m ; n decreases as Z mul cx r -1 / 2 . Even in 



6 



Takeuchi & Ida 



10" 



N 



£10"' 



10" 



(a) T =150K "" 


0.01 


M^=0.01 M„. ir , 






0.03. 


r Mh^-osm.,. 






" M disk =u. I M^_ a 




0.1 • 




L oc r 



0.1 



10 



100 



f(AU) 



10" 



N 



10"' 



(b) 7" =300K 


0.01 




0'.03 


mL=0.03 M 






M dil=°' 1 M sun 






L oe r 



0.1 



10 



100 



r(AU) 



Fig. 7. — Minimum dust abundances Z m i n for various disk mod- 
els, (a) Cold disks (Tq = 150K). (6) Hot disks (T = 300K). Red 
solid and blue dashed lines represent models with shallow and steep 
density profiles, S 9 oc r~ 1 and S g oc r~ 1,5 , respectively. 

the hot gas disks (To = 300K) with a mass as small as 
-Mdisk = 10~ 2 Af Q , the required dust abundances for SGI 
are less than 0.1. If the temperature at the disk midplane 
is as cold as To S3 150K, as expected in passive disks that 
are heated only by stellar radiation (Chiang & Goldreich 
1997; Tanaka et al. 2005), SGI operates for Z min < 0.05 
at 10AU of disks with M disk = 10 _2 M Q . 

4. DtSCUSStON 

Analysis in this paper assumes that turbulence is in- 
duced in the dust layer. If the dust abundance at the 
midplane / m jd is larger than unity (Z > (Cog?)) 1 / 2 ), tur- 
bulence weakens with increasing Z as a oc Z~ 2 (Equa- 
tion flSJ). The growth timescale of SGI rapidly decreases 
as T grow oc Z~ 4 (Equation ©). A slight increase in 
Z makes SGI viable. Yll also estimated the minimum 
dust abundance Z n required for SGI (Equation (69) of 
Yll), assuming that the turbulent diffusion parameter 
of a dense dust layer (/mid ^ 1) was a r] ~T s f] (Equation 
(68) of Yll), independent of Z. This estimate for a n 
adopted by Yll corresponds to the maximum value of 
a(Z) at Z ss [Cesf]) 1 / 2 in our estimate (Equation (|5j|). 
If a had a constant value a v for all the condition 
TU-ow < Tdrift would be 



Z > Z v ~ r) 



5/6 2/3 
'9 



(14) 



where we used / m id = Zh g /hd and hd/h g ~ rj 1 ^ 2 . Note 
that this Equation ((14"]) differs from Equation (69) of Yll, 
because we use Equation f|7|) . in which Thrift oc ,/midi 
while Yll assu mes Thrift oc /^ id (see discussion following 



Equation (IA6|) below). As discussed in T12, accumula- 
tion of the dust increases the inertia of the dust layer 
and decelerates the inward drift of the dust, resulting in 
weaker turbulence. Considering a oc Z~ 2 , the condition 



L grow 



< Tdrift becomes 



Z> Z„ 



2/SqX/S 



(15) 



as shown in Equation (|10p . We ignore numerical co- 
efficients such as C c ff and C n . These two conditions 
(Equations fli"4"l) and ([15]) ) coincide when / m ;d = 1 (i.e. 



^ - ~mm ~ V 1/2 ) or when Q g = Q g ,o(f]) = i] 
(i.e., on the dashed line in Figures EUSJ . If the disk's 
Q g is greater than Q ffl o i then the latter condition (Equa- 
tion (]T5]) ) permits smaller dust abundances for SGI by 
a factor of Z min /Z n = (Qg/Qgfi)^ 1 ^ 3 ■ For example, at 
1AU on a disk with Mdisk = O.O3M , S g oc r _1 , and 



,-i/2 



T = 300K, Equation JI5J gives Z min = 0.11, which is 
smaller by a factor of 2 than the value Z v = 0.23 pre- 
dicted by Equation (TT4]) . Considering the Z dependence 
of a reduces the required minimum value Z m i n for SGI 
as compared with that derived by Yl 1 . 

5. SUMMARY 

We analyze the condition for SGI occurrence in the 
dust layer. The growth timescale of the instability must 
be shorter than the radial drift timescale of the dust. 
The growth timescale decreases as turbulence weakens. 
The necessary condition for SGI is obtained by consider- 
ing turbulence induced in the dust layer. Using the tur- 
bulence strength estimated from the energy supply rate 
from the accreting dust, the minimum dust abundances 
for SGI, Z m i n , are derived as a function of the Toomre 
Q g parameter of the gas disk and the deviation fraction rj 
of the gas velocity from the Keplerian value. If the dust 
particles are small and their stopping time is less than the 
Keplerian time (T s < 1), Z m ; n is independent of T s . For 
disks with Q g - 10 2 and 77 - 10~ 2 , SGI occurs if Z > 0.1. 
The required Z decreases with decreasing Q g and 77, be- 
coming as small as 0.01 for Q g ~ 10 and 77 ~ 10~ 3 . Such 
an increase in Z from the solar abundance is expected 
to occur through several processes, including the radial 
drift of the dust and dispersal of the gas from the disk 
(Youdin & Shu 2002; Takeuchi & Lin 2002; Takeuchi et 
al. 2005). Therefore, SGI provides a possible route for 
small particles to directly collapse to planetesimals in an 
initially laminar disk. If the gas disk were globally tur- 
bulent via MRI, for example, the required Z would be 
larger than Z m i n estimated in this paper. 
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Fig. 8. — Comparison of approximate expressions with numer ical result s for (a) turbulent diffusion parameter a and (b) dust radial 
drift velocity i^r- Solid lines represent approximate expressions l)A4j l and ||A6| |. and dashed lines are numerically calculated by T12. (a) 
Diffusion parameter a is plotted for T a = 10 — 3 ,10 -2 , and 10 — 1 . (b) Radial drift velocity r is plotted for T s = 10 — 3 . 



APPENDIX 

APPROXIMATE EXPRESSIONS OF DIFFUSION PARAMETER AND RADIAL DRIFT RATE 

In this Appendix, we evaluate an approximate expression of the turbulent diffusion parameter a for small particle 
limit (T s <C 1) according to T12. We consider turbulence excited via KH or streaming instabilities of the dust layer. 
The strength of the turbulence is controlled by the dust accretion rate. The dust particles accrete toward the star 
either due to the gas drag on individual particles (Nakagawa et al. 1986) or collective drag exerted on the entire dust 
layer (Weidenschilling 2003; Youdin & Chiang 2004). Individual drag dominates if the dust-to-gas ratio in the dust 
layer is less than unity (/mid ^ l)i while collective drag is important if / m id > 1. 

First, we consider individual drag (for / m id <j 1). Using the formula derived by Nakagawa et al. (1986; or Equation 
(9) of T12, hereafter Equation (T9)), the accretion velocity of small dust particles (T s < 1) is 



Vd,r 



-2T s i]v K 



(Al) 



To satisfy angular momentum conservation, the gas moves outward with the velocity v g , r — —Vd, r Pd/ ' Pg, where pd (p g ) 
is the dust (gas) density. The effective gravities (including the pressure gradient force) acting on the dust and gas are 
gd = — rClfr and g g = —(1 — 27y)rf7^, respectively. The net work done by these effective gravities on a unit surface area 
of the dust layer is calculated as AZ?drag ~ / (pdgdVd,r + Pg9g v g,r)dz <~ ri v^l-K.T s T,d (Equation (T12)). The energy 
dissipation due to turbulence is estimated as AEturb ~ Si ay erUe ddy /Teddy, where £i aye r is the surface density of the dust 
layer, and u c ddy and r e ddy are the typical values for the velocity and turnover time of the largest eddies, respectively. 
Adopting the so-called "a prescription", we evaluate w c ddy ~ \foLC g and r c ddy ~ Or (Cuzzi et al. 2001). For /, m d < 1, 
the column density of the dust layer is dominated by the gas, and thus £i aycr ~ ^ghd/h g . Here, hd/h g is the ratio of 
the scale heights of the dust layer and gas disk, and is given by hd/h g yJa/T s (Equation (T4)). Then, the energy 
dissipation rate is A_E t urb ~ \J a 3 /T s h 2 g Vt^Yi g . Equating the energy supply and dissipation rates of turbulence, with 
an efficiency parameter for the energy supply, AE'turb ~ C e g A^dragj gives 



a 



(C eS fjZ) 2 ^T s 



(A2) 



The collective drag force in the ^-direction on a unit surface area is estimated by the plate drag approximation for 
/mid ^ 1 as Pq z r*> p g vdv g fi / 'dz ~ — p g vr]VK / hd, where v is the turbulent viscosity, and the velocity shear in the vertical 
direction dv g fijdz is estimated as "qv-^/hd- This drag causes accretion of the dust layer and an outward motion of the 
upper gas layer. The dust layer loses energy ^l^rPg z and the upper gas layer gains energy (1— rj)QKrPg z , where the work 
done by gas pressure is taken into account. The net energy liberation rate is AE^s 



(Equation (T18)), where we used v 



hjT s Q,K and hd 



dissipation rate is Ai? tU rb 
Ai?turb ~ C e ff AEvis gives 



J layerW cddy 



/^ddy ~ ah 2 Sl^E d , 



g J mid , 

where Si aye r 



-T]fl K rP ez ~ T] 2 v^n K T s T, d /f m id 
(Equations (T16) and (T6)). The energy 
~ is dominated by the dust. Equating 

(A3) 



where we use / m id ~ Z^/T s /a (Equation (T6)). Expressions (|A2|) and (|A3|) are the same as Eq uati on (T28), except for 
numerical factors and the adoption of a slightly steeper power-law index (S = 1 in Equation (|A3|) instead of 6 = 0.94 
in Equation (T28)) to simplify analytical expressions. Connecting these expressions gives 



(A4) 



Takeuchi & Ida 



where the numerical fac tors C% — 1.0 and C2 = 1.6 are set to fit the numerical result of T12. In Figure [SJa), the 
approximate expression (IA4[) is compared to the numerical result of T12. The approximation overestimates a around 



Z = \/C c ffri and underestimates it for large Z. However, the error is less than 2 for 10 < Z < 10 , and the overall 
behavior appears to be acceptable. 

The midplane dust-to-gas ratio / m id = Zh g /hd ~ Z^/T s ja is, from Equations (|A2I) and (|A3I) . / m id ~ [Z 2 /(Cegfj)] 1 / 3 
for / m id ^ 1, and / m ;d ~ Z 2 / (C e sfj) for / m id ^ 1- Connecting these two expressions gives an approximate formula for 
/mid as 

fuM « ( ttIWT + ( T^r-\ ■ (A5) 



CiCcsfjJ \C 2 C e sfi 

From the above Equation (|A5|) . it is seen that / m ;d ~ 1 when Z s» (Ccfffj) 1 / 2 . 

The drift velocity of the dust in the limit of / m id <C 1 and T s <C 1 is given in Equation (jAll) . In the limit of / m id 3> 1, 
Vd.r ~ rPg z / '(vK^d) ~ V v kT s / f m id (Equation (T51)). Connecting these limits, we estimate 

2T s 77u K 

/mid + 1 

which gives the dust drift timescale given in Equation ([7]). This approximate expression is compared with the numerical 
calculation by T12 in Figure Hth), which shows that the above expression provides a good approximation. Note that 
v d,r oc f~l d for / m ;d 1. Equation (2.11) of Nakagawa et al. (1986) shows that, due to individual drag, the radial 
velocity of the dust particles at the midplane is proportional to f~? d - However, the average velocity of the dust layer, 
/ pd,Vd,rdz /Ed, is approximately proportional to /~jd, provided that the dust density profile pd{z) is Gaussian. The 
drift velocity due to collective drag is also proportional to /~;d, as shown in Equation (|A6j) . Thus, for / m id > 1, the 
radial drift time derived from the particle velocity at the midplane (Td r ;ft oc f 2 lid ), which is adopted in Yll, SC11, 
and MKI12, overestimates the average drift time (Td r ;ft oc /mid); resulting in a less strict condition for T grow < Tdrift 
than reality. 
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